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Abstract: An improved algebraic reconstruction technique(ART) with tunable diode 
laser absorption spectroscopy(TDLAS) is presented in this paper for determining 
two-dimensional distribution of H2O concentration and temperature in a simulated 
combustion flame. This work aims to simulate the reconstruction of spectroscopic 
measurements by a multi-view parallel-beam scanning geometry and analyze the 
effects of projection rays on reconstruction accuracy. It finally proved that 
reconstruction quality dramatically increases with increasing number of projection 
rays until they are more than 180 for 20X20 grid, after that point, the number of 
projection rays has little impact on reconstruction accuracy. It is clear that the 
temperature reconstruction results are more accurate than the water vapor 
concentration by traditional concentration calculation method, the study in this 
article also proposed an innovative way to remove the error of concentration 
reconstruction and improve the reconstruction quality greatly, the capability of this 
new method is evaluated by appropriate assessment parameters. By using this new 
approach, the concentration reconstruction accuracy is not only badly improved, but 
also a suitable parallel-beam arrangement is put forward for sake of high 
reconstruction accuracy and simplicity of experimental validation. At last, a bimodal 
structure of combustion region is assumed to demonstrate the robustness and 
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universality of the proposed method. Numerical investigation indicates that the 
proposed TDLAS tomographic algorithm is capable for accurate temperature and 
concentration profiles detection, this feasible formula for reconstruction research is 
expected to resolve several key issues in practical combustion devices. 

Key words: algebraic reconstruction technique, two-dimensional distribution, 
projection rays, accurate reconstruction 


PACS: 42.62.Fi, 02.70.-c 


1. Introduction 

Over the past decades, TDLAS has become an indispensable method for making 
sensitive and fast in situ measurements in gaseous detection, such as species 
concentration, pressure, temperature, and velocity[1,2,3]. Traditional absorption 
spectroscopy is a line-of-sight(LOS) technique, which can only reflects averaged 
value along the laser path without providing any spatially resolved information[4,5,6], 
while the temperature and concentration among combustion area are usually 
non-uniform, thus, it is necessary to accurately obtain the spatial distribution of fluid 
field parameters. Recently, TDLAS data is combined with computer tomography(CT), 
herein referred to as tunable diode laser absorption tomography(TDLAT), to image 
two-dimensional distribution of combustion properties of scramjet wind tunnel. 
TDLAT can be roughly divided into two categories in terms of reconstruction 
algorithm: transformation method, the representative algorithm is filtered back 
projection(FBP) algorithm[7,8,23], which needs complete projection data to realize 
physical reconstruction, and the other one is iterative method, the main example is 
algebraic reconstruction technique(ART) algorithm[9,10,11,12]. Bryner et al. 
employed filtered back projection algorithm to reconstruct the peak of the water vapor 
absorption transition centered at approximately 7164.9cm’' of flat burner and indoor 
environment based on experimental measurements[13]. Li Fei Et al. designed a set of 
TDLAT system with 6 parallel rays which were arranged in different position, they 
reconstructed the temperature and concentration distribution of scramjet combustor 


exit using ART[14]. L. Ma used a tunable laser called “fiber Fabry-Perot tunable filter 


laser(FFP-TFL)” to scan over the 1333-1377nm range at ~200Hz in the 2D 
reconstruction of temperature and HO concentration distributions over a Hencken 
burner flame[15]. Although the FBP algorithm can obtain highly accurate 
reconstructions of objective area, especially in medical imaging, a large number of 
closely sampled angular views and projection rays may increase the difficulty of 
experimental design, therefore, FBP reconstruction study based on laboratory 
verification is rarely reported in the literature. Unlike FBP, ART tomography 
algorithm has been successively used for combustion diagnosis because of its strong 
resistance to noise and availability with only limited projection data[12,16]. Even 
though many researchers have performed relevant work for imaging, it is hard to 
evaluate the preciseness of reconstruction results, especially the concentration 
distribution results, it is usually not so satisfactory when compared with temperature 
reconstruction by using traditional method[23]. 

In this paper, we divided the combustion area into 20x20 grids, and theoretically 
discussed the reconstruction of temperature and concentration with different 
projection rays based on a modified adaptive algebraic reconstruction 
technique(MAART)[17]. More importantly, we proposed a novel way to eliminate 
the concentration reconstruction error, which greatly improved the accuracy of the 
results when compared with the traditional method, furthermore, an optimal spatial 
structure of test system include parallel-beam projection rays and grid partition 
principle was suggested later. At the same time, asymmetric temperature and 
concentration distribution models were created to analyze the robustness and 
universality of the proposed algorithm. By the way of this article, the temperature and 
concentration can be simultaneously and precisely generated whatever the combustion 


flow field is like in the future. 


2. Theoretical foundation 
Generally the TDLAT can be divided into two steps, which is absorption 
spectroscopy theory on the one hand and tomographic reconstruction algorithm, on 


the other. 


2.1 Absorption Spectroscopy Principle 
According to the Beer-Lambert law, when collimated light at a certain frequency 
penetrates a gas medium of path length, L, the relationship between the transmitted 
and incident intensity of radiation is: 
1,(V) =V] Xa ()SIT 16,4) a) 
Where I, and I, are the transmitted and incident laser intensity respectively, p is 


total static pressure[atm], X» is the mole fraction of the absorbing species, S,(T) is 


abs 


the transition line strength[cm™~ atm] for a single transition i, and ¢,is the 
normalization line-shape function[cm], namely, f ¢,dv =1, usually ø, approximates 


the convolution of two dominant broadening mechanisms, i.e., Doppler and 
collisional broadening, representative for typical combustion environment. 


The integrated absorbance of the transition i can be inferred from the above as: 


9S, [T (x)]dx (2) 
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The absorption coefficient can be defined as flows from equation(2) for a uniform 
flow field whose concentration and temperature are constant: 

a= PX Si (T) (3) 
The line strength S of a single transition i, which is temperature-dependent, can be 


described in terms of the formula below at a reference temperature T): 
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Where h [J.s] is Plank’s constant, c [cm.s7] the speed of light, k [J.K] the 
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Boltzmann’s constant, Q(T)the partition function of the absorbing species[18], v, 
[cm"] the line-center frequency, and E' [cm"] the lower state energy of the transition. 
2.2 Tomography algorithm 

Before the tomography algorithm study, we discrete the region of interest(ROI) into 
N(N equals to nxn) grids, and the flame parameters such as pressure, temperature and 
gas concentration are assumed to be constant in each grid, then the integrated 
absorbance A of transition iin equation(2) for each laser beam is calculated by the 


discrete form: 
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Where I and J are respectively total beam and grid number, L, the path length of the 


i-th laser beam passing through the grid j, which can be automatically calculated from 
a fast algebraic reconstruction algorithm based on improved projection coefficient 


computation[19], @,, the absorption coefficient of the i-th laser beam within the j-th 


vj 
grid. In general, equation(5) can be rewritten as matrix expression: 


A=La (6) 
Where L is the IXJ matrix of weights L,, the column vector @ = {@,1, Œ, 25- Œ y Y 


and A={4 1, 42 AF. 


Instead of the matrix theory, here we compute œ by using a modified adaptive 
algebraic reconstruction technique(MAART)[17], application of this ART method 
solves linear equations(6) in iterative equation as follows: 

ak =a" -A la".L-A] (i=1,2,...,1) (7) 
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Where k represents the iteration time in the ART procedure, 4, is the relaxation 


A, = BX (j=1,2,...,J) (8) 


parameter, and £ generally is a constant( # <0.5)[17] during the calculation, here we 
make it 0.25. The iterative process is terminated when the absorption coefficient in 
equation(3) becomes less than 1X 10° between two consecutive iterations. 
Considering the non-negativity of concentration and temperature value, restrictions 


are performed during the iterative process: 


[®(@)], -| i j (9) 


In order to reduce the mutation of reconstruction results, the reconstructed value is 
to be restricted within a reasonable range for specific application in practice, here a 
kind of smoothness regularization is employed[21,22], for grid (i,j), the absorption 


coefficient a is: 
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The a. p after k-th iteration is not only relative with (k-1)-th grid information, but 


also dependent on its nearest eight neighbors, and we keep the constant 6 0.001[17]. 


By performing the equations(7). (8). (9), the tomographic reconstruction result of 
acan be obtained subsequently, the temperature of grid j would be approximately 


determined by the following formula if we select a pair of appropriate transitions(v, 


and v,) with different temperature dependence[20]: 
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The traditional way to infer concentration is that X, can be simultaneously obtained 
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from equation(3) at the atmospheric pressure with T, in hand, but this algorithm may 


generate great reconstruction error, we propose a novel method instead of this way to 


solve the problem, more details will be given in the following text. 


3. Methods 
3.1 Absorption line selection 

H20 is an attractive target species for combustion diagnosis due to the relative 
abundance of water vapor as a combustion product and a strong near-infrared 
absorption spectrum. In this work two transitions 7153.722cm"&7153.748cm™ and 
7154.354cm'', have been selected to infer temperature by means of two-line scheme 
because of their suitable line-strength and free of significant interference from nearby 
transitions[9,23], these two overlapped transitions 7153.722cm™ and 7153.748cm™ 
can be treated as one transition when we deal with their line strength, integrated 
absorbance, center frequency . The dependence of the two line strength on 
temperature can be obtained from HITRAN 2008 or be calculated by equation(4), as 
shown in Fig.1(a), then the line strength ratio of the two absorption transitions and 


relative sensitivity to temperature are derived as shown in Fig.1(b). 
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Figure.1(a). Line strength for selected transitions 
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Figure.1(b). Ratio and sensitivity 


3.2 Numerical investigation for temperature reconstruction 

In the simulation, we assumed the ROI is 20x20 cm°, then the ROI was divided 
into 20X20 grids, namely, the length for both X and Y axis directions of each grid was 
1cm, which was a relatively high spatial resolution to understand combustion flow 
field information. Parallel-beam retained 20 light rays at each projection angle, and 
these rays rotated around the ROI for generating projection data in different views. This 


parallel laser-detector scanning geometry is shown as Fig.2. 
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Figure.2. Parallel beam collection geometry for multi-view with a fixed location 


A Gaussian model was designed for simulating the temperature and concentration 
distribution of the combustor exit of flat flame furnace, and we can arbitrarily specify 
a possible distribution using this Gaussian function. In this paper an interest 
temperature range(500K-1400K) was taken for example, similarly, water vapor 
concentration i.e. volume fraction varied from 0.15 to 0.25 for testing, and the total 
static pressure was always latm, which is the case for many practical applications. 
First of all, projection data of the numerical phantom, which is exactly equivalent to 
the integrated absorbance A [26], can be obtained directly from equation(5). 
Secondly, we changed the number of projection angle which is scattered between 0 
degree and 180 degree to bring out the effect of projection rays on the overall 
reconstruction quality by the algorithm mentioned in this article, in the end, we 
specifically set up the projection rays varying from 120 to 360 with 60 increment and 
keep the parallel rays 20 invariably at each projection angle. Take one of beam 
arrangements for example, the reconstructed profile was shown in Fig.3, here we just 


displayed the temperature result regardless of concentration, which would be studied 


later. In order to analyze the impact of projection rays on reconstruction accuracy 
clearly, rather than other influence factor, we did not add measurement noise to 


projection data. 
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(a). The original temperature distribution (b). The reconstructed temperature distribution 
Figure.3. Temperature image reconstruction for Gaussian model 
by ART with totally 300 projection rays scattered between 0 degree and 180 degree 


As can be seen from Fig.3, very good agreement was achieved between the 
reconstruction image and the original phantom, considering the other beam 
arrangements mentioned in the above, the maximum relative error between the 
tomographic images and the phantoms was utilized to quantitatively evaluate the 
accuracy of the reconstructions, this criteria was mathematically expressed as[25]: 


lo, -r 


o = Onin 
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e = Max, (12) 


Where o,and r, represent the true temperature value and reconstructed one for grid 
j respectively, N represents the total grid number in reconstruction area, and o,,,. 
represents maximum value of original pixel, similarly, o,,, the minimum, e reflects 


the maximum reconstruction difference of a certain pixel compared to original 
phantom. Another criterion for evaluating reconstruction quality is the average error, 
shown in equation(13)[25],which inflects average difference of each pixel between 


reconstructed image and original one. 
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In terms of the different parallel-beam projection rays scanning, temperature 


reconstruction error is plotted in Fig.4(a) by calculating equation(12) and (13) . 
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Figure.4. Reconstruction error versus projection rays number 
From the graph above we can obviously conclude some typical trend is: the error 


decreases with increasing projection data, reconstruction accuracy is significantly 
dependent on the projection rays when they are less than 180, after this point, 


projection rays have little impact on e and aver e, which indicates that there is no 


need to increase the projection rays as many as possible in the TDLAT experiment, for 
20x20 grid partition, in consideration of laboratory verification and experimental 
condition, 180 or near this magnitude of projection rays arrangement is accurate 
enough to reconstruct the temperature distribution. 
3.3 Modified algorithm for concentration reconstruction 

Just as described in most of literature, once the temperature of each grid is known, 
the line-strength S can be calculated, and then the local concentration follows from 
equation(3) as: 

X,=a,/ pS(T,) (14) 

In terms of the water vapor concentration reconstruction involved in this paper, we 

analyzed the results just like temperature mentioned in the above, the concentration 


image was mapped in Fig.5 under the same conditions as temperature. 
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(a). The original concentration distribution (b). The reconstructed concentration distribution 
Figure.5. Concentration image reconstruction for Gaussian model 
by ART with totally 300 projection rays scattered between 0 degree and 180 degree 
Fig.5 shows the water vapor concentration reconstruction result under a specific 


parallel-beam layout, we can also plot the trend curve of the maximum relative error 
e and the average error avere versus different projection rays by equation(12) and 
(13), which was displayed in Fig.4(b). 

It turns out that the concentration can be easily reconstructed as long as the 
temperature has been imaged correctly. However, it is evident that the concentration 
reconstruction results are usually worse than the temperature reconstruction either 
from Fig.4 or Fig.5, especially at the edge of ROI, reconstructed value of 
concentration generates relatively lager error than the center region, this phenomenon 
can be explained by the literature[23]. Many researchers have been focused on 
temperature measurement for a long time, but they pay little attention on the 
correction of concentration reconstruction. 

Currently, an unique and standard way of solving the concentration distribution is 
through directly iterative calculation instead of indirect derivation by temperature 
reconstruction results. In practical measurements, with T in hand, the line strength 
S(T) for each grid can be figured out easily, then we defined another projection 
weights L', called normalized L, which is depicted by the multiplication of 
original weights L at each grid and the corresponding line-strength, shown as: 

L'=L-S(T) (15) 

Since the pressure of the combustion region remains latm at the beginning of our 
study, integrated absorbance A for every transition of a certain laser beam will be 
converted to equation(16), which would be called absorption formula for transition i. 

A=X,,,°L' (16) 

According to the previous text, there would be 2I absorption formulas because 

the total projection laser beam number is I, now we combine all these 27 equation 


together. 
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Where X, represents the concentration of grid j, L,,.,, the normalized path 


length of the i-thlaser beam within the j-th grid with regard to the first absorption 


transition vl, and Lyg 


for another, A„; and A,, represent integrated 
absorbance of absorption vl andv2 respectively for i-th laser beam. Equation(17) 
is a linear equation array with J unknowns and 2I equations, solution of the 
equations is deeply restricted because of the increased I equations when compared 
with equation(5) , which can facilitate us to obtain a more accurate outcome than 
before. As we all know from the matrix theory, there must be at least x equations to 
solve x unknowns if we want to get the unique and precise solution, therefore, it 
provides a reference guideline for us to select the appropriate laser projection beam 
arrangement in TDLAT, in other words, projection rays should approximate the half of 
grid number if we want to get accurate results with small number of beam scanning as 
possible as we can. 


In order to reconstruct the concentration distribution by ART from equation(17), we 


treated X, as iteration variable and the final convergence value was exactly 


concentration. In combustion diagnosis, there could be a wide variety of flow field 
model, a relatively large range(0.15 to 0.25) of water vapor concentration has been 
tested in the above, now another concentration layout varying from 0.15 to 0.1591, 
which is an almost uniform distribution, has been supposed by changing the 
parameters of Gaussian function to assess the capability and accuracy of the proposed 
method, for the sake of simplicity, temperature value was still the reconstructed 
results listed previously. Fig.6 shows the reconstruction for 20X20 grid by traditional 
method(TM) and the modified method(MM) on the condition that the combustion 


zone is scanned by 300 projection rays . 


0.16 0.16 
0.155 0.155 
0.15 0.15 
0.145 ash 


(a). H2O concentration of reconstruction by TM (b). Interpolation operation result 


(c). Original concentration distribution (d). Interpolation operation result 
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Figure.6. Reconstruction image of assumed concentration distribution 


In order to obtain a high spatial resolution image, cubic spline interpolation 
algorithm was used to realize a 200x200 grid distribution, which were displayed in 
Fig.6.(b). (d), (f). Itis evident from Fig.6 that traditional method generates much 
larger mistake than modified one, especially at the edge, obvious fluctuation exists 
whether the original reconstruction or the interpolated image shown in Fig.6(a) and 
Fig.6(b). While the proposed calculation scheme implements extremely well, the 
reconstructed results is almost the same as original image shown in Fig.6(e) and 
Fig.6(f). To quantitatively compare the reconstruction accuracy between TM and MM, 
we calculated the maximum relative error e and the average error avere 
respectively by using equation(12) and (13) again, a list of the calculation results is 


shown in table 1 for comparison. 


Table 1. Error analysis for 300 parallel-beam scanning 


Reconstruction maximum relative average error 
methods error €e aver e 
Traditional 
method(TM) 1.186 0.1643 
Modified 


method(MM) 0.0673 0.0054 


It can be seen that a significant improvement of the tomographic accuracy of 
concentration reconstruction can be obtained by using the proposed method. To 
further validate the modified concentration calculation algorithm, we employed this 
modified way to conduct the reconstruction whose ROI is the same as the distribution 


depicted by Fig.5(a) again. Figure.8 shows the concentration reconstruction result. 
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(a). Original concentration distribution (b). Reconstructed concentration distribution by MM 
Figure.7. Concentration reconstruction result by modified concentration calculation algorithm 
with totally 300 projection rays scattered between 0 degree and 180 degree 


In some complicated combustion environment, the very specific case of the 
centered symmetrical Gaussian distributions of temperature and mole fraction 
discussed above should be expanded to include a range of possible cases, we 
considered a bimodal structure to validate the universality and robustness of the 


proposed method, reconstruction results are displayed in Fig. 8. 
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(e). Temperature reconstruction residuals (f). Concentration reconstruction residuals 
Figure.8. Reconstruction of asymmetry structure for a bimodal distribution 


It is obvious from Fig.5, Fig.7 and Fig.8 that the modified concentration 
calculation algorithm can not only reconstruct different kinds of concentration 
distribution easily but also provides a more accurate reconstruction image than the 
traditional method, especially at the edge of ROI, irregular fluctuation of 
reconstructed phantom is greatly eliminated. In general, the modified method for 
reconstructing concentration distribution is basically a universal approach in practical 


combustion diagnosis. 


4. Conclusion 

In summary, a comprehensive analysis for high-quality and simultaneous 
tomographic imaging of temperature and gas concentration was implemented in this 
paper. We numerically demonstrated a spatial resolution with totally containing up 
20X20 grid points for an improved ART algorithm. By quantitatively calculating the 
reconstruction error of different parallel-beam scanning arrangements, it proved that 
reconstruction accuracy was greatly dependent on projection rays only when they 


were less than a certain value, additionally, the concentration reconstruction results 


are always worse than the temperature reconstruction by traditional algorithm. 
Therefore, a modified concentration calculation method was introduced later, which 
aimed to directly reconstruct concentration instead of deriving from temperature. This 
proposed method facilitates the high-quality concentration reconstruction without 
adding extra projection rays. The conclusions presented in the above indicate that 
projection rays should be approximately equal to half of the grid number or near this 
magnitude if we simultaneously consider the reconstruction accuracy and simplicity in 
engineering application. At the end of this article, we applied this proposed method to a more 
complicated flow field model for demonstrating the robustness and universality of the method. 
The study in this work validates an available and feasible approach for precisely 
reconstructing temperature and concentration distribution, and we believe that it certainly can 


offer a good theoretical foundation for experimental verification in the future. 
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